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Skymaps measured with imaging atmospheric Cherenkov telescopes (lACTs) represent the real 
source distribution convolved with the point spread function of the observing instrument. Cur¬ 
rent lACTs have an angular resolution in the order of 0.1° which is rather large for the study of 
morphological structures and for comparing the morphology in y-rays to measurements in other 
wavelengths where the instruments have better angular resolutions. 

Serendipitously it is possible to approximate the underlying true source distribution by applying a 
deconvolution algorithm to the observed sky map, thus effectively improving the instruments an¬ 
gular resolution. From the multitude of existing deconvolution algorithms several are already used 
in astronomy, but in the special case of y-ray astronomy most of these algorithms are challenged 
due to the high noise level within the measured data. 

One promising algorithm for the application to y-ray data is the Maximum Entropy Algorithm. 
The advantages of this algorithm are the possibility to take a priori knowledge into account and 
that it is an independent approach to previous work, e.g., [1] who applied the Richardson Lucy 
Algorithm to y-ray skymaps. 

An implementation of the Maximum Entropy Algorithm is provided in the MemSysS software 
package [2]. As this algorithm is very sensitive to various input parameters it is essential to un¬ 
derstand their influences. We present a study of the influences of these parameters in order to in¬ 
vestigate the applicability of the Maximum Entropy Algorithm for the deconvolution of skymaps 
in y-ray astronomy. 
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1. Introduction 

Nature does not reveal its mysteries easily. This principle is demonstrated by e.g. studying 
morphological structures. Due to the limited angular resolution of detectors, images represent 
always the source distribution convolved with the point spread function (PSF). In the case of very- 
high-energy (VHE; >100 GeV) y-ray images (skymaps) of imaging atmospheric Cherenkov tele¬ 
scopes (lACTs) the size of the PSF is in the order of 0.1°, which is rather large, especially compared 
to other wavelengths, and thus causes a few challenges. 

Hence even a perfect point source exhibits an extent - however exactly the question whether 
or not a source is a point source is often of crucial interest, cf. for instance [3]. 

A large PSF also limits the possibilities to study structures within the source, especially in the 
order of or smaller than the PSF. A paramount example for this is the differentiation whether a 
supernova remnant (SNR) exhibits a shell type morphology or not. 

In addition multiwavelength studies are essential in y-ray astronomy to obtain further inside 
into the processes taking place in astrophysical sources. This study is biased by the fact that the 
PSF of y-ray experiments is often orders of magnitude larger compared to other wavebands. For 
instance the PSF of the X-ray observatory XMM-Newton is two orders of magnitude better than the 
one of H.E.S.S. To overcome this problem the traditional ansatz is to smooth data from the superior 
instrument with the PSF from the other one, obviously loosing information by doing this. 

One way to diminish these problems is the application of adequate deconvolution algorithms 
to VHE y-ray skymaps to approximately derive the underlying true source distribution, thus effec¬ 
tively improving the instruments angular resolution. From the multitude of existing deconvolution 
algorithms several are already used in astronomy, but in the special case of y-ray astronomy most 
of these algorithms are challenged due to the high noise level within the measured data. In previous 
works algorithms like the Richardson-Lucy Algorithm [1] have successfully been applied to VHE 
y-ray data. In this work we will report a study about the applicability of the Maximum Entropy 
Algorithm to VHE y-ray data. 

The main object of the study, presented in this work, is the investigation of parameter influ¬ 
ences on the deconvolution and finding an optimal parameter set. This study is based on simula¬ 
tions. In Section 2 the simulation of y-ray excess maps is outlined. In the subsequent Section 3 the 
Maximum Entropy Deconvolution, as used in this work, will be described, followed by an overview 
over our main results in Section 4. To conclude we will give a brief summery and outlook in Sec¬ 
tion 5. 

2. Simulation of y-ray skymaps 

To simulate a y-ray excess map realistic input parameters e.g. for the camera acceptance, the 
background rate and the PSF have to be chosen. In our case templates were obtained from the 
H.E.S.S. Model++ analyses [4] of the supernovae remnant (SNR) RX J1713.7—3946 and the 
Crab nebula for extended and for point sources, respectively. The main background source in y-ray 
astronomy with Cherenkov telescopes are events from Cosmic rays which are evenly distributed 
over the field of view. Therefore a homogeneous background distribution is assumed. The num¬ 
ber of background events is taken from the analysis mentioned above as well as the number of 
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excess events. The latter one is scalable according to the simulated source significance. As source 
morphology different source type can be chosen, e.g. point like. 




Figure 1: Schematic view of the simulation process (not to scale). Left panel: Flat background (dashed 
red) with a source (dashed blue) on top of it. In solid colours the background and source modulated with 
the 7 -ray acceptance. Central panel: The black bars represent the simulated 7 -ray and background events. 
Right panel: Simulated events corrected for the 7 -ray acceptance. 

The main steps of the simulation process are illustrated in Figure 1 and are given below. 

1) The background and the source distribution is convolved with the 7 -ray acceptance. The 
dashed lines in Figure 1 {left panel) show the unconvolved background {red) and source 
{blue) distributions. The solid areas give the convolved ones. 

2) Random events are generated according to the distribution obtained in the first step. This is 
illustrated in Figure 1 {central panel), the black lines illustrate the simulated events. 

3) The map obtained in the second step is corrected for the acceptance. As illustrated in Fig¬ 
ure 1 {right panel) the source in the centre is only slightly affected by the correction, but the 
background distribution resembles a flat distribution with larger fluctuations at the edges of 
the skymap. 

4) In the final step the above assumed background level is subtracted. 

In Figure 2 examples of different source morphologies and simulated skymaps are given. The 
models of extended sources are adjusted to resemble the 7 -ray emission of RX J1713.7—3946. 

3. Maximum Entropy Deconvolution 

The Maximum Entropy Method [5] is widespread in a surprisingly large variety of fields. Its 
applications range from studying the impact of the climate change on coffee supply chains [ 6 ] 
to calculations of drug absorption rates [7]. In astronomy it is used in nearly all wavelengths 
regimes: from radio astronomy (cf., e.g. [ 8 ]) to optical images from the Hubble Space Observatory 
(especially in the period before the refurbishment of the detector system, see, e.g., [9, 10]) up to 
the 7 -ray observatory COMPTEL [11]. 

Its basic principle is that the most realistic solution is the one which minimises the amount 
of information while remaining compatible with the data. This minimisation of information is 


3 






















































Application of Maximum Entropy Deconvolution to y-ray Skymaps 


Susanne Raab 





11'’15"'00® 11''10'"00" 11'’15"'00® H^IOW 11''15"’00' 11''10"’00" 

Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 


Figure 2: Top row: Examples of source models, which can be used for the simulation. From left to right: 
A point source, a Gaussian, a ring and a modified ring shape, all fitted to resemble the 7 -ray emission of 
RX J1713.7—3946. Bottom row: Simulations according to source shapes displayed in the top row. Besides 
the point source on the left hand side, the simulations are smoothed with a Gaussian (a = 0.05°) for better 
visual perception. 


equivalent to maximising the entropy. From this theory a method was invented to reconstruct 
images from incomplete and noisy data, cf. [12]. In this work Maximum Entropy Deconvolution is 
used in the implementation provided by the MemSysS software package [2]. 

The deconvolution procedure has several input parameters, which are listed in Table 1. 


Parameter 

Description 

Configurations 

maxiter 

Number of iterations 

1 , 2 , 4,..., 1024 

nscales 

Number of hidden images, each progressively blurred, to be 
combined to form the visual image 

7,8 

icfwidth 

Width of the intrinsic correlation function 

0, 1,2,...,8 

default 

Value of default distribution 

0.001,0.01,0.1, 1 

aim 

Value of the stopping criterion 

0.001,0.01,0.1, 1 

noise 

Standard deviation of the Gaussian noise 

0.0001 

n s amp1e 

Number of samples taken to calculate standard deviation etc. 

3 

bayes 

Stopping criterion; “1” refers to classic maximum entropy 
stopping criterion 

1 

entropy 

Type of entropy; “1” refers to standard entropy 

1 


Table 1: List of parameters used for studying parameter infiuences on the deconvolution. Eor an in-depth 
description of their definition see [2]. The configurations denote the tested parameter combinations. 


To obtain deep inside into the algorithm and to achieve the best possible results it is essential 
to understand their influence on the deconvolution. We have performed a detailed study of the 
influence of the input parameters. Therefore we performed a multidimensional scan with 3168 
parameter sets (for the applied steps see Table 1). To obtain statistically reliable results this scan 
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was repeated ten times for different simulations. Examples of deconvolved skymaps can be seen 
in Figure 3. As the deconvolution result proved to be robust against varying values of icf width 
only the results of the remaining four parameters are shown in Section 4. 




05 '^ 36 '^ 00 ® 05 ^ 34 ^ 00 ® 05 ^ 32 '^ 00 ® 

Right Ascension (J2000) 



Figure 3: Examples of deconvolved skymaps of a point source with a significance of 50 sigma. From left to 
right: The values for the parameter maxi ter are 1, 8, 64 and 1024, respectively. 


4. Analysis & Results 


To judge which parameter set yields the best deconvolution results, we follow the example 
of [1] and take the relative error (cf. Equation 4.1) as a measurement for the goodness of the 
deconvolution. 

H [^deconvolution (t j) ~ Amodel (t 7*)] 


RE = 


hj 




a Amodel (t 7) 
ij 


(4.1) 


For the deconvolutions of a point source there is a second method to appraise the results: a 
Gaussian is fitted to the deconvolved skymaps, the width of the Gaussian is a measure for the 
angular resolution and can therefore be used as a criterion to find the optimal parameter set. 
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Figure 4: Results of a systematic scan of deconvolutions for point sources of significances, colour coded 
from 10 sigma (blue) to 100 sigma (green). From top to bottom, the parameters maxiter, nscales, 
default and aim are given. In the left hand column the relative error RE normalised to the REsim of 
the undeconvolved simulations is given, on the right hand side the width cr of a Gaussian fitted to the 
deconvolved skymaps, normalised to the width (Tsim of the undeconvolved simulations. The grey areas 
depict the standard deviations of REsim and (Tsim- 
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We evaluated the results of parameter scans for point sources with significances between 10 
and 100 sigma. The best parameter set, e.g. that yielding the smallest values for the relative error 
RE and the Gaussian widths, was taken to illustrate the influence of the parameters in Figure 4. 
All parameters were kept fixed at their best values except the one shown in the respective image. 
On the left hand side the relative errors RE are given, normalised to the relative error REsim of the 
undeconvolved simulations. On the right hand side the same for the width of a Gaussian fitted to the 
deconvolved skymaps. The errorbars each give the standard deviation from the ten deconvolutions 
using the same parameter set. 

As seen on Figure 4 {top panel) the parameter maxiter has the largest influence and the 
width of the Gaussian can be improved for values of maxiter > 64 by up to a factor of 1.95 
compared to the one of the undeconvolved simulations. Larger values might be able to slowly 
improve the results further, but are accompanied by an approximately linear increase of computing 
time. The other parameters, illustrated in Figure 4 {secondary to bottom panel), are harder to grasp. 
Due to the larger error bars no definite statement can be made, but one might perceive some trends. 
The parameter aim seems to favour smaller values, just like the parameter default. 

What can clearly be seen on all graphs, is that deconvolution works generally better for sources 
with a higher significance and can improve the angular resolution by a factor up to 1.95. Further¬ 
more the optimal parameters tend to be the same for all tested significances. 

5. Conclusion & Outlook 

The detailed study, presented here, has shown, that the Maximum Entropy Deconvolution is a 
feasible choice to improve the angular resolution of VHE y-ray skymaps which yield a rather bad 
signal-to-noise ratio compared to other wavelengths. 

The potential improvement in the angular resolution is compatible with the results of a similar 
study with an independent algorithm [1]. In addition it can be asserted, that it is important to 
understand the influences of various parameters in order to derive a good set of parameters before 
applying the deconvolution algorithm to real data. 

The next logical step will be the systematic study of parameter influences for extended sources 
of different types and the application to real data. 
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